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Resonant inelastic x-ray scattering (RIXS) tuned for the Cu L edge is a possible tool to detect 
charge excitations in cuprate superconductors. We theoretically investigate the possibility for ob¬ 
serving a collective charge excitation by the RIXS. The RIXS process via the intermediate state 
inevitably makes the spectral weight of charge excitation stronger in electron doping than in hole 
doping. Electron-hole asymmetry also appears in the dynamical charge structure factor, show¬ 
ing a new enhanced small-momentum low-energy mode in electron doping. These facts indicate a 
possibility of detecting the new charge mode by RIXS in electron-doped systems. 

PACS numbers: 78.70.Ck, 78.20.Bh, 74.72.-h 


I. INTRODUCTION 

The cuprate superconductor is classified as a doped- 
Mott insulator. The spin and charge degrees of freedom 
in the cuprates, thus, exhibit their characteristic dynam¬ 
ics. Recently, resonant inelastic x-ray scattering (RIXS) 
experiments tuned for the Cu L edge have provided a lot 
of new insights about their dynamics. [l| In hole-doped 
cuprates, paramagnon excitations near the magnetic zone 
boundary are less damped [2| unlike the theoretical pre¬ 
diction of the t-t'-t"-J model considered to be a canon¬ 
ical model of cuprates. Q In an electron-doped cuprate 
Nd 2 -ccCea,Cu 04 (NCCO), the energy of magnetic exci¬ 
tations identified by RIXS as well as inelastic neutron 
scattering increases with carrier concentration Hi in 
contrast with hole-doped systems. The hardening of the 
magnetic mode is explained by taking into account a 
three-site correlated hopping process ignored in the t- 
J-type model. [6[ 

The charge ordering in the pseudogap state of hole- 
doped cuprates is a recent hot topic, to which RIXS 
has contributed significantly. 0, [![ Even for the electron- 
doped NCCO, the presence of charge ordering has been 
indicated by resonant x-ray scattering. [9| A charge¬ 
ordering signal is enhanced below a temperature where 
antiferromagnetic (AFM) fluctuations are developed. 
The RIXS experiments for NCCO (Refs. 4 and 5) have 
reported a mode that disperses higher in energy than the 
magnetic one. Whereas the mode has been interpreted 
as a signature of broken symmetry based on its temper¬ 
ature dependence [5| similar to the charge ordering, it 
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has been identified as a charge excitation Q because of 
the agreement with incoherent intraband charge excita¬ 
tion observed by Cu K-edge RIXS. [Mini If the mode 
is related to charge ordering, a similar one should be ob¬ 
served in L-edge RIXS for hole-doped cuprates, but such 
a mode has not been reported so far. Therefore, clarify¬ 
ing the difference in charge dynamics between hole- and 
electron-doped cuprates seen by L-edge RIXS is crucial 
for putting forward our understanding of cuprates. In 
Cu K-edge RIXS, such a difference has been attributed 
to the difference of the screening process. [T2| In con¬ 
trast to K-edge RIXS, there is no theoretical confirma¬ 
tion about the difference of L-edge RIXS between hole- 
and electron-doped systems. 

In this paper, we theoretically investigate electron-hole 
asymmetry of charge excitations seen by L-edge RIXS, 
focusing on two possible contributions: One is from the 
RIXS process where charge fluctuation is affected in the 
intermediate state with the Cu 2p core hole, and the other 
is the difference of dynamical charge structure factor it¬ 
self. For the former contribution, we start with a two- 
band model of the Cu02 plane and calculate the RIXS 
spectrum by using the exact diagonalization technique. 
The charge spectral weight in the hole-doped system is 
suppressed in RIXS through the intermediate state where 
hole carriers are influenced by the core hole. On the other 
hand, the weight in the electron-doped system is insen¬ 
sitive to the intermediate state because of no electron 
carrier on the core-hole site. This behavior is also seen 
in the single-band Hubbard model. These facts natu¬ 
rally explain visibility of charge dynamics in electron- 
doped cuprates. For the latter contribution, we per¬ 
form a large-scale density-matrix renormalization-group 
(DMRG) calculation of dynamical charge structure factor 
in the Hubbard model with next-nearest-neighbor hop- 
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ping. We find a new enhanced small-momentum low- 
energy charge excitation in the electron-doped system. 
The energy of the excitation is lower than that of spin 
excitation. This excitation mode is different from the pre¬ 
viously confirmed incoherent intraband excitations whose 
energy is higher than spin excitation. Combining the two 
contributions, we predict that the L-edge RIXS will de¬ 
tect the new charge collective excitation in the electron- 
doped cuprates. 

This paper is organized as follows. The RIXS spec¬ 
tra of the two-band model and the single-band Hubbard 
model are shown in Sec. mi for both the electron and 
hole dopings. In Sec. Hill we calculate the dynamical 
charge structure factor of the single-band Hubbard model 
obtained by dynamical DMRG and discuss a new en¬ 
hanced small-momentum low-energy charge excitation in 
the electron-doped system. Finally, a summary is given 
in Sec. IIVI 


II. RIXS SPECTRUM 


We firstly consider a two-band model of the Cu02 
plane that contains a Cu 3 (ia^ 2_^2 orbital and a O 2p 
Wannier orbital with x‘^ — symmetry, jl^ which is 
equivalent to a standard three-band model of cuprates 
from which a nonbonding band is removed. Since we 
focus on Cu L-edge RIXS associated with the Mx 2 _y 2 
orbital, the use of the two-band model can be justified. 
The Hamiltonian reads 

Hsd = 2Tpd y] Tij (dl^(j)ja + H.C. j 

ija 

+ ^ X X ( 1 ) 

ia i 

where the operator dia- is the annihilation operator for 
a ?idj. 2 _y 2 hole with spin a at site i, = d\^dicr^ (pia 
represents the annihilation operator of the symmetric O 
2p orbital, Tij = ^ (3^^ with the position 

vector Ri and /3k = [sm‘^{kx/2) + sin^(/cy/2)]“^/^, A is 
the level separation between 3 da^ 2_^2 and 2p orbitals, and 
Ud is the Coulomb repulsion of the Cu 3 da, 2_^2 orbital. 
We take Tpd = 1 eV, A = 3 eV, and Ud = S eV, which 
are typical values appropriate for cuprates. 

In RIXS, tuning the polarization of incident and outgo¬ 
ing photons, we can separate excitation with the change 
of total spin by one (AS = 1) and excitation with no 
change of total spin {AS = 0) [l3-[lq (see Appendix [A|) . 
The two excitations can be defined as 


with S{ = and 

1 


oi _ -iq-Rz jf 


I 


- Hi 


Eq + iT 


dicj, (4) 


where |0) and |/) represent the ground state and final 
state with energy Eq and Ef^ respectively; j is the to¬ 
tal angular momentum of Cu 2p with either j = 1/2 or 
j = 3/2; R^ is the position vector at site /; ooi is the 
incident photon energy; T is the relaxation time of the 
intermediate state; and = H^d + with 

Uc and Sj being the Cu 2p-3d Coulomb interaction and 
energy level of Cu 2p, respectively. Here, we assume the 
presence of a core hole at site 1. We note that, if we take a 
fast-collision approximation assuming T is very large, [l| 
and become the standard charge operator Aq and 
the 2 ; component of spin operator of a 3 dx 2 _y 2 hole, 
respectively. 

In order to calculate Eqs. m and m, we use the 
Lanczos-type exact diagonalization technique on a ^/l0 x 
VTO-unit-cell cluster under periodic boundary condi¬ 
tions. We set Uc = 9 eV and T = 0.5 eV. The incident- 
photon energy Si is tuned to the energy of the Cu 2p ab¬ 
sorption peak that is calculated separately (not shown). 
Figure [11(b) shows foi* fhe 

carrier concentration x = 2/10 = 0.2 in hole doping. 0 
The low-energy excitation below 1 eV is dominated by 
the AS = 1 excitation, while the AS = 0 excitation is 
higher in energy. These behaviors are consistent with 
calculated RIXS of a single-band Hubbard model. [H 
The results of the x = 0.2 electron doping are shown in 
Fig. Ha). The distribution of and is similar 

to the hole-doped case, but the spectral weight of 
is larger in the electron-doped case. This means that the 
charge excitation is more visible in electron doping than 
in hole doping. 

The electron-hole asymmetry of partly comes 

from the RIXS process where charge fluctuation is af¬ 
fected in the intermediate state with Cu 2p and Cu 3d 
Coulomb interaction. In order to clarify this, we calcu¬ 
late integrated spectral weight up to 2 eV 

for the AS = 0 {AS = 1) excitation and take the ratio 
Jr = We plot the ratio at q = (37r/5,7r/5) 

as a function of Uc in Fig. |2](a). [TsI In hole doping, 
Jr decreases with decreasing Uc showing a minimum at 
Uc ^ A = 3 eV, while Jr is less dependent on Uc in 
electron doping. Smaller charge spectral weight in hole 
doping shown in Fig. [1] (14 = 9 eV) is related to the 
decrease of Jr as explained below. 

We also performed the exact diagonalization calcula¬ 
tion of for the single-band Hubbard model 

given by 


(A:^) = X K/l Hi |0) f <5 {Aco -Ef + Eo) (2) 

/ 

(Aco) = X I (/I |0) 4 - % + Ho) (3) 

/ 


Htu = X 4aCi+<5<T + y; ni^nn, (5) 

idcr i 

where is the creation operator of an electron with 
spin a at site i, number operator riia- = i E S 
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FIG. 1. (Color online) RIXS spectra on the ten-unit-cell clns- 
ter of the two-band model O- (a) Electron doping and (b) 
hole doping with x = 2/10 = 0.2. Black and red colors repre¬ 
sent AaS = 1 and 0 excitations, respectively. The solid curves 
are obtained by performing a Lorentzian broadening with a 
width of 0.2 eV on the delta functions denoted by vertical 
bars. Parameters are = 8, A = 3, Tpd = 1, f/c = 9, 
r = 0.5 in units of eV. Incident photon energies for RIXS are 
set to the edges of absorption spectra. 





FIG. 2. (a) The core-hole potential Uc dependence of 

spectral-weight ratio /r = at q = (37r/5,7r/5) 

for the ten-unit-cell cluster of the two-band model, where Is 
are obtained by integrating the weight up to Aa; = 2 eV 
in Fig. [1] Filled and open diamonds represent the x = 0.2 
electron- and hole-doped cases, respectively, (b) The same as 
(a) but for the 18-site single-band Hubbard with U = lOt and 
r = t at q = (7r/3,7r/3). The weights are integrated up to 
Acc; = 6t. Filled and open diamonds represent x = 2/18 = 
0.11 electron-doped and hole-doped cases, respectively. 


represents the four nearest-neighbor sites around site i, 
and t and U are the nearest-neighbor hopping and on-site 
Coulomb interaction, respectively. In calculating L-edge 
RIXS, we introduced an attractive Coulomb interaction 
Ucc between the core hole and the valence electron. 

Figure [2fb) shows Ucc dependence of the ratio Ir 
defined by the integrated weight up to Acj = 6t at 
q = (7r/3,7r/3). With decreasing Ucc down to Ucc ^ t/. 


Ir in hole doping decreases, which is similar to Fig.[2]^a). 
Since there is no electron-hole asymmetry in the single¬ 
band Hubbard model, the asymmetric Ir dependence is 
purely caused by the effect of Ucc working at the interme¬ 
diate state. In electron doping, the incident photon only 
kicks the site where there is no electron carrier, i.e., the 
undoped site, since an electron-doped site has no empty 
state to accommodate a core electron. In this case, the 
intermediate state of RIXS does not create additional 
charge modulation, leading to small Ucc dependence on 
Ir. In hole doping, we have to take into account both 
undoped and hole-doped sites. When the undoped site is 
excited by the photon, double occupation of electrons is 
induced with energy cost U — 2I/cc- On the other hand, 
the hole-doped site gives the energy cost —Ucc- In cal¬ 
culating RIXS, we tune the incident photon to the site 
with lower energy cost. When Ucc — ^ 0, the undoped 

site is excited as is the case of electron doping, leading to 
charge fluctuations similar to the electron doping. With 
decreasing Ucc^ the hole-doped site starts to participate 
in the intermediate state. This eventually suppresses Ir 
through destructive contribution to the hole-doped sites. 
The same type of discussion can be applied to the two- 
band model shown in Fig.[2fa). The fact that the AS = 0 
excitation is larger in the electron-doped system than in 
the hole-doped system is thus caused by the RIXS pro¬ 
cess related to the core-hole Coulomb interaction. We 
note that the rapid increase of Ir for Ucc < U is caused 
by the decrease of spin excitation. 


Next let us consider the effect of the asymmetric elec¬ 
tronic state between hole- and electron-doped cuprates 
on RIXS. Such an asymmetric state is introduced by 
electron hopping beyond the nearest-neighbor sites. We 
note that such a long-range hopping is induced by di¬ 
rect oxygen 2p-2p hopping. Including the next-nearest- 
neighbor hopping term given by -t'a 4a^i+S'^ be- 
ing S' the four next-nearest-neighbor sites, into the single¬ 
band Hubbard model ([5]), we calculate the RIXS spectra 
for a ^/Ts X VTs periodic lattice with two electrons and 
two holes {x = 0.11) as shown in Fig. [3l [l3] Here, we 
take t' = —0.25t. Comparing Fig. [Sfa) and Fig. [^b), 
we find that the AS = 0 excitation at q = (7r/3,7r/3) 
is enhanced for electron doping. at q = (7r/3,7r/3) 
in Fig. Ha) is 0.68, which is larger than the case of 
t' = 0 with Ir = 0.55. This increase can be attributed 
to the asymmetric electronic state. We thus conclude 
that both the RIXS process and the asymmetric elec¬ 
tronic state contribute to enhancing charge excitations 
in the electron-doped system seen by RIXS, being con¬ 
sistent with the fact that charge excitations have been 
reported in NCCO (Ref. 4) but not yet in hole-doped 
cuprates though similar charge excitations should exist. 
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FIG. 3. (Color online) RIXS spectra on the 18-site single¬ 
band Hubbard cluster with U = lOt and t' = — 0.25t, Ucc = 
15t, and F = t. (a) Electron doping and (b) hole doping with 
X = 2/18 = 0.11. Black and red colors represent AS = 1 and 
0 excitations, respectively. The solid curves are obtained by 
performing a Lorentzian broadening with a width of 0.2t on 
the delta functions denoted by vertical bars. Incident photon 
energies for RIXS are set to the edges of absorption spectra. 


III. DYNAMICAL CHARGE STRUCTURE 
FACTOR 



FIG. 4. (Color online) Dynamical structure factors in a 6 x 6 
cylindrical Hubbard cluster with t' = — 0.25t and U = St. 
Spin structure factor S'(q,ct;) at q = (tt/T, tt/S) for (a) elec¬ 
tron doping and (b) hole doping. Charge structure factor 
iV(q,ct;) at various q for (c) electron doping and (d) hole dop¬ 
ing. Black, red, green, purple, and blue lines represent spectra 
with X = 0 (half filling), 2/36, 4/36, 6/36, and 8/36, respec¬ 
tively. A Gaussian broadening width of 0.2t is used for the 
spectral weights. 


For making a prediction to RIXS, it is important 
to clarify characteristics of low-energy charge excitation 
coming from asymmetric electronic states. We thus cal¬ 
culate dynamical charge structure factor A(q, cj) for a 
6x6 single-band Hubbard cluster with cylindrical geom¬ 
etry, where the x direction is of open boundary condition 
while the y direction is of periodic boundary condition, 
by using the ^namical DMRG method. We use U = 8t, 
t' = —0.25t. 0 The numerical method and accuracy are 
detailed in Appendix [0 

In order to confirm doping dependence of spin dynam¬ 
ics, [il we show in Figs.|ll^a) andlU^b) the dynamical spin 
structure factor 5'(q, cj) for electron and hole doping, re¬ 
spectively, at q = (tt/T, 7r/3). Whereas the peak position 
slightly shifts downward with hole doping, the position 
moves to higher energy with electron dopin^as reported 
by the quantum Monte Carlo calculations, p 

Figure m^c) exhibits the doping and momentum depen¬ 
dence of A(q, cj) for electron doping. At q = (tt/T, 7r/3) 
that is the smallest q defined in the cluster, the low- 
energy spectral weights below uj = 1.6t show an inter¬ 
esting doping dependence: The weight increases with in¬ 
creasing X from x = 0.06tox = 0.11 but it decreases with 
further increase of x followed by an upper shift of spectral 
weight. More interestingly, the main spectral weight is lo¬ 
cated in the energy region below spin excitation centered 
around uj = t. This is a characteristic behavior near half 
filling in the dop ed Mott insulator with strong AFM spin 
correlation. This also means that this charge mode 
is different from incoherent intraband charge excitations 
that have been observed in RIXS experiment for NCCO 


as a broad and largely dispersive structure above mag¬ 
netic excitation. Q In this sense, the charge mode shown 
in Fig. m^c) has not yet been identified in experiment. 
Since its spectral weights may overlap with those of spin 
excitations in RIXS experiments, as is expected from our 
calculations, improving resolution and/or a polarization 
analysis for outgoing photons is necessary to identify this 
charge mode in RIXS experiment. 

Such a low-energy charge mode is not pronounced in 
hole doping as shown in Fig. SJ^d). This contrasting be¬ 
havior may be related to the strength of spin correlation 
in the spin background. In fact, as is the case of the t- 
J-type model, the AFM spin correlation is large in 
electron doping than in hole doping near half filling as is 
seen by the static spin structure factor (see Appendix ICj) . 
The large AFM spin correlation in the electron-doped 
system can make a strong spin-charge coupling that in¬ 
duces a collective mode of charge assisted by spin, [2T| 
while such an effect is weak in hole doping. 

The presence of the low-energy charge mode assisted 
by ^in has been discussed by using t-J-type mod¬ 
els. |3. [21| - [^ However, its strength depending on carriers 
has not been discussed by using a realistic parameter set 
for cuprates. Furthermore, there are only a few works on 
dynamical charge structure factor of the Hubbard-type 
model. [MH We thus emphasize that the present re¬ 
sults clarify for the first time the low-energy charge mode 
of electron-doped cuprates by using a realistic Hubbard- 
type model. 
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IV. SUMMARY 

In summary, we have theoretically investigated 
electron-hole asymmetry of charge excitations seen by L- 
edge RIXS and clarified two contributions to the asym¬ 
metry. One is related to the RIXS process, in which the 
charge spectral weight in the hole-doped system is sup¬ 
pressed through the intermediate state where hole carri¬ 
ers are influenced by the core hole, while the weight in 
the electron-doped system is insensitive to the presence 
of the intermediate state. Another contribution is due to 
the asymmetric electronic state inducing a large charge 
spectral weight in electron doping. These facts natu¬ 
rally explain visible charge excitations in electron-doped 
NCCO, but not in hole doping. To characterize low- 
energy charge excitation, we have performed a large-scale 
DMRG calculation of dynamical charge structure factor 
in the Hubbard model with the next-nearest-neighbor 
hopping. We have found an enhanced small-momentum 
low-energy charge excitation in the electron-doped sys¬ 
tem. Based on these findings, we predict that the L-edge 
RIXS will detect the new charge collective excitation in 
the electron-doped cuprates. 
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Appendix A: RIXS formula 

The spectral weight of Gu L-edge RIXS with momen¬ 
tum transfer q and energy transfer Acj is given by [iHi3 

/q (Au;) = ^ |(/| Tq |0)|2 5(Au; - Ej + Eo), (Al) 
/ 

with Tq = (at-^q + wherc |0) and |/) rep- 

resent the ground state and final state with energy Eq 
and Ef, respectively, j is the total angular momentum 
of Gu 2p with either j = 1/2 or 3/2. The operator 
induces the change of total spin by one {AS = 1 ), 


but does not charge the total spin {AS = 0 ). The 
operators are given in the main text. The coefficients 


a- 


N 


and in Tq are determined by the polarization 


vector of the incident photon ei and the outgoing pho- 

= 2(eo 


ton Co 


a 




eQef)/15 and 


= 2i{eo X ei)Vl5. 


Appendix B: Dynamical DMRG calculation of 
dynamical structure factors 


The dynamical charge and spin structure factors, 
V(q, cj) and 5'(q, cj), are defined as 

«(q,.) = -ilm (01 |0XB1) 

SM = -ilm(0| 10) (B2) 

where and S^ are the Eourier component of charge 
operator ni = ni^ ^4 and that of the ^ component of 
spin operator S[ = ( 77,4 — 77 , 4 )/ 2 , respectively, with the 
number operator nia- at site I and spin <j; T is the system 
size; H is the Hamiltonian of a given system; and 7 is a 
small positive number. 

We calculate Eqs. m and (Irnli for the Hubbard 
model with the next-nearest-neighbor hopping t' by dy¬ 
namical density-matrix renormalization-group (DMRG) 
method. We use a a 6 x 6 cluster with cylindrical geom¬ 
etry where the x direction is of open boundary condition 
while the y direction is of periodic boundary condition. 
In the cluster {Lx = 6 and Ly = 6 ), the y component of 
momentum q is determined by using a standard transla¬ 
tional symmetry, i.e., qy = Tiy'KjLy {uy = 1 , 2 ,-•• ,T^), 
but the X component is given by qx = nx'K/{L J- 1) 
{ux = 1 , 2 , • • • ^Lx) because of the open boundary condi¬ 
tion. Defining lx and ly as the x and y component of site 
/, respectively, we can write 


^ y {Lx + 1)U^ 


■ sin(^ 3 ^/ 3 ^) e ^ Til 




(B3) 

(B4) 


In dynamical DMRG, we use three target states for 
a given uj: Eor ^(q,^^) ([Eq. IBlj . (i) |0), (ii) N'qlO), 
and (hi) {uj — H Eq 27 ) |0). The target state 

(iii) is evaluated by using a kernel-polynomial expansion 
method. [25| In our kernel-polynomial expansion method, 
the Lorentzian broadening 7 is replaced by a Gaussian 
broadening with half width at half maximum 0 . 2 t, t being 
the nearest-neighbor hopping amplitude of the Hubbard 
model. In our numerical calculations, we divide the en¬ 
ergy interval [ 0 , 2 t] by 100 -mesh points and have targeted 
all of the points at once. 

To perform DMRG, we construct a snakelike one¬ 
dimensional chain, and use the truncation number m = 








6 



FIG. 5. (Color online) The dependence of N{c[,uj) for a give 
parameter set on the truncation number m in the dynamical 
DMRG method. A 6 x 6 cylindrical cluster with t' — — 0.25t, 
on-site Coulomb interaction f/ = 8t, and x = 4/36 for electron 
doping is used, q = (tt/T, tt/S). A Gaussian broadening width 
of 0.2t is used for the spectral weights. 

2000, and resulting truncation error is less than 2 x 10“^. 
In order to check the accuracy of our results, we examine 
m dependence of A^(q, cj). Figure [5] shows the results for 
the carrier concentration x = 4/36 at q = (tt/T, 7 r/3). It 
is clear that the difference between m = 1500 and 2000 
is small, though the convergence is not perfect. Since 
it takes 16.5 h by using 3500 nodes in the K computer, 
RIKEN Advanced Institute for Computational Science, 
to obtain a single curve of A/'(q, uj) with m = 2000 by our 
dynamical DMRG code, it is practically difficult to in¬ 
crease m further to perform systematic calculations of the 
doping and momentum dependence of dynamical struc¬ 
ture factors. 


Appendix C: DMRG calculation of static spin 
structure factor 

In the t-t'-t"-J model, it is known that antiferromag¬ 
netic spin correlation remains strong in the electron- 
doped side of the phase diagram. [2^ This is caused by 
the same sublattice hopping t' affecting the spin back¬ 
ground through the stability of the Neel-type configura¬ 
tion by the motion of electron carriers. The antiferromag- 
netically ordered spins can couple to the charge degree 
of freedom, leading to the low-energy collective charge 
excitation in the t-J-type model. In order to check 
whether our t-t'-U Hubbard model examined by dynam¬ 
ical DMRG has a similar effect on the spin background, 
we investigate the static spin structure factor defined by 

S{q) = (0| |0), (Cl) 

which can measure the strength of antiferromagnetic spin 
correlation. 

Figures [6] and [7| show S{q) for the electron- and hole- 
doped systems, respectively. It is clear that near half¬ 
filling for the electron-doped side [Figs. [6l[a) and [6l[b)] 
the antiferromagnetic spin correlation is dominating in 
the spin background as evidenced from the large value of 
S{q) near q = (7r,7r). This is in contrast with the hole 
doped case in Fig. [71 We consider that the contrasting 
behavior of spin correlation contributes to the different 
behavior of low-energy charge dynamics discussed in the 
main text. 

Another contrasting behavior appears in the over¬ 
doped region at x=0.22, where an incommensurate spin 
correlation developed in the hole-doped case (Fig. [Tfd)) 
but not in the electron-doped case (Fig. |6l[d)). The dif¬ 
ference is probably related to the Fermi-surface topology 
that favors a nesting in the hole-doped side. 
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FIG. 7. (Color online) The static spin structure factor *S'(q) for the hole-doped 6x6 cylindrical Hubbard cluster with t' — —0.25t 
and on-site Coulomb interaction U = St. (a) x=0.06, (b) x=0.11, (c) x=0.17, and (d) x=0.22. The height of the bars represents 
the magnitude of S{q) at a given momentum. 






































